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ABSTRACT 

We measure X-ray emission from the outskirts of the cluster of galaxies PKS 0745-191 with 
Suzaku, determining radial profiles of density, temperature, entropy, gas fraction, and mass. 
These measurements extend beyond the virial radius for the first time, providing new infor- 
mation about cluster assembly and the diffuse intracluster medium out to ~ 1.5r2oo (?'200 — 
1.7 Mpc ~ 15'). The temperature is found to decrease by roughly 70percent from 0.3— lr2oo- 
We also see a flattening of the entropy profile near the virial radius and consider the impli- 
cations this has for the assumption of hydrostatic equilibrium when deriving mass estimates. 
We place these observations in the context of simulations and analytical models to develop a 
better understanding of non-gravitational physics in the outskirts of the cluster. 

Key words: galaxies: clusters: individual: PKS 0745-191 - X-rays: galaxies: clusters - galax- 
ies: clusters: general 



1 INTRODUCTION 

The outskirts of galaxy clusters present an opportunity to study the 
formation of large scale structure as it happens. Beyond the core 
where non-gravitational processes such as cooling flows and feed- 
back from active galactic nuclei can dominate activity, clusters are 
expec ted to be mor e relaxed and to follow self-similar scaling rela- 
tions l l Kai ser" 19861. The virial radius can be thought of as a border 
between regions of equilibration and infall, so merger activity from 
accreting material may play an important role in the outer dynam- 
ics. 

As the largest virialized systems in the universe, clusters of 
galaxies can be useful in constraining cosmological models. The 
position of clusters at the high end of the mass spectrum and their 
evolution from the initial density perturbations makes them sensi- 
tive probes of the scale of these fluctuations, erg, and the density 
of matter in the universe, Qm ( White et al. 1993; Eke et al. 1996). 
Additionally, the evolution with redshift of cluster properties such 
as the gas mass fraction can be use d to study other p ortions of the 
cosmic energy density budget (e.g., lAUen et al .12003) . For a review 
of the cosmological importance of galaxy clusters, see j2005l) . 

To understand cluster properties, it is important that the re- 
lationships between observables and derived quantities are well- 
calibrated. A common example is that temperature and den- 
sity profiles can be measured from X-ray spectra and used to 
determine a cluster's mass under the assumptions of spherical 
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symmetry and hydrostatic equililbrium between the gas pres- 
sure and gravitational potential. Though numerical simulations 
predict a decline in temperature towards the virial ra dius (e.g., 
lEvrard et al.ll 19961 : iFrenk et al Jl 19991 : iLoken et al J[2002l) . observa- 
tions at smaller radii have produced inconsi stent results. Some anal- 
yses have found decli n ing pr o files (e.g., 



Markevitch et al. 1998 



De Grandi & Molendil |2002|: IVikhlinin et al.. ,2005. : .Pratt et al 



2007h . while others have seen a scatt ering of slopes consistent with 
flat or even increasing profiles (e.g., Ilrwin et al.lll999l : lzhang et al.l 
2004; Arnaud et al. 2005). These studies typically do not measure 
the temperature profiles much further out than half the virial radius, 
even with Chandra and XMM. Different methods of estimating the 
virial radius, using either scaling relations or direct determination 
from the mass profile, make comparisons difficult, but very few 
measurements have b een made of cluster gas temperatures out t o 
the virial radius (e.g., ISolovveva et"ai]|2007l; iReiprich et all2008l) . 
and to our knowledge none have been reported beyond it. 

Until recently. X-ray detectors have been unable to measure 
the temperature of the intracluster medium (ICM) out to the virial 
radius due to its low surface brightness in the outskirts relative 
to background noise. Typicall y a mass profile such as that of 
iNavarro, Frenk & White! ( 1 19971 NFW) is used to extrapolate mass 
estimates to distances greater than those observed. One aim of our 
work is to measure properties of the ICM out to the virial radius 
in order to check the reliability of such extrapolations. For our pur- 
poses, we equate the virial radius with r2oo, the radius within which 
the mean total density is 200 times the critical density of the uni- 
verse at the redshift of the cluster. 

The low orbit of Suzaku places it within Earth's mag- 
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netopause, giving it a significantly lower and more stable 
particle background compared to Chandra and XMM-Newton. 
iReiprich et all dlOOSl) have recently leveraged this ability for low 
surface brightness observations in the outskirts of A2204, measur- 
ing the temperature nearly to r2oo- Other cluster observations with 
Suzaku have demonstrated its capacity for t emperature and abun- 
dance measurements (e.g., ISato et aUbOOTl) . providing an outline 
for some of the methods used here. 

In this work, we present Suzaku observations of PKS 0745- 
191, a relaxed, c ool core cluster tha t is the brightest i n X-rays be- 
yond z = 0.1 jpabian et alj|l985l : lEdge etaLlll99(]| : lAUen et all 
Il996l) . Previous measurements have found a mean gas temperature 
in the range of 6 .4 — 8.5 keV, depending on the models used and 
regions studied (Chen, Ikebe & Bohringeij l2003l and references 
therein). [Allen et al. (1996.) find good agreement between the cen- 
tral mass from X-ray s and that de termined from a strongly lensed 
arc, but the results of IChen et al.l ( 120031) disagree, finding a factor 
of 2 smaller mass from XMM data (see also [Snowden et al. 20Qi|, 
for a reanalysis of the XMM data). We describe our observations in 
the next section, followed by details of the spectral analysis in Sec- 
tion[3l and the resulting profiles in Section|4l In Section|5]we place 
our findings in the context of other cluster studies, and we conclude 
the paper in Section|6] For distance scales, we use the cosmological 
parameters Ho = 70kms~^ Mpc~^ and fim = 1 — ^^a = 0.3, 
giving an angular scale of 113 kpc per arcminute at the cluster's 
redshift, z — 0.1028. The plotted and quoted error ranges are la 
statistical uncertainties except where otherwise stated. 



2 OBSERVATIONS AND DATA REDUCTION 

Suzaku observations of PKS 0745-191 were taken between 2007 
May 11-14 in five separate fields of roughly 32 ks each. Details 
of these pointings are listed in Table [T] The central pointing is 
aimed toward the peak of the cluster emission and the others are 
positioned adjacently, with 4' overlap at each chip edge, as in 
Fig. [T| We use only data from the X-ray Imaging Spectrometer 
(XIsT fKovama et al JlioOTi) . We exclude the back-illuminated de- 
tector, XISl, despite its higher sensitivity at low energies, because 
it also has a significantly higher particle background level than the 
two front-illuminated sensors, XISO and XIS3, which have similar 
responses. There is an offset between the measured temperatures of 
the BI and FI chips, and we discuss this contribution to our system- 
atic uncertainties in Section 1431 We simultaneously analyze data 
from the two FI sensors to effectively double the exposure time 
while keeping the noise low. Observations were taken in normal 
clocking and editing modes with spaced-row charge injection on, 
and the data have been processed with the energy scale calibration 
of v2.1.6.16. The data preparation described below was carried out 
using XSELECT 2.4.1 and FTOOLS version 6.5.1 from HEASARC, 
with instrumental parameters from CALDB updated 2008 Septem- 
ber 5. Updates to the calibration parameters released during the 
preparation of this paper do not affect the results beyond the statis- 
tical uncertainties. 

We performed the standard screening of events files to re- 
move time intervals during satellite manoeuvres, telemetry satu- 
ration, and passage through the South Atlantic Anomaly, as well as 
elevation angles less than 5° and 20° from the nighttime and day- 
time Earth. The light curve for the remaining time intervals is sta- 
ble, showing no sign of flaring or excess particle background. We 
exclude the corners of the detectors where ''''Fe calibration sources 
lie, and obvious point sources seen in Suzaku or XMM images were 



Table 1. Observational parameters of the five pointings 



Obs. ID 


Position 


Exposure (ks) 


RA 


Dec (J2000) 


802062010 


Center 


32.0 


116.8852 


-19.2901 


802062020 


NW 


32.2 


116.6543 


-19.2063 


802062030 


NE 


30.8 


116.9737 


-19.0727 


802062040 


SE 


32.9 


117.1155 


-19.3739 


802062050 


sw 


33.4 


116.7966 


-19.5079 




Figure 1. Exposure-coirected mosaic of PKS 0745-191, smoothed with a 
Gaussian of radius 1'. Background subtraction and vignetting correction 
have been omitted, and calibration regions and point sources have not yet 
been excised. Ring radii are 2.5', 6', 9.5', 13.5', 18.5', and 24'. 



excised with circular regions of radius 2.5' to remove more than 
99 per cent of the flux spread by the PSF. Detector response ma- 
trices and effectiv e area functions are constructed with xisrmfgen 
and xissimarfgen dlshisaki et al.ll2007l) . respectively. The latter tool 
generates the spectral response for a given sky brightness, account- 
ing for the known issue of contamination on the optical blocking 
filter. Uncertainties in this contamination should not affect our re- 
sults significantly, as we do not consider energies below 1 keV. Vi- 
gnetting effects will be different for cluster and background emis- 
sion components, so we use a uniform surface brightness out to 
a radius of 20' when modeling the background emission and in- 
put a /3-model su rface brightness profile using the parameters of 
IChen et al. I( l2003l) for the auxiliary response file applied to the clus- 
ter emission model. We find that modifying the parameters of the 
surface brightness profile used to determine the "effective area" 
with xissimarfgen does not significantly impact the output. Thus, 
we do not expect an unresolved central cusp or uncertainties from 
the outward extrapolation of the surface brightness profile to influ- 
ence our results. 



© 0000 RAS, MNRAS 000, 000-000 



PKS 0745-191 at the virial radius 3 



2.1 Background Subtraction and Modelling 

Accounting for background emission is critical when observing re- 
gions of low surface brightness, and it consists of multiple com- 
ponents. We subtract the non X-ray back ground (NXB) of charged 
particles and gamma rays with xisnxbgen l lTawa et alj2008h . which 
models the trend of these events from night Earth data, weighted by 
the magnetic cutoff rigidity. The X-ray background consists of solar 
wind charge exchange, thermal emission mainly from the hot local 
bubble and the Galactic halo, and the cosmic background (CXB) 
due to unresolved point sources. Charge exchange contaminates the 
low energy spectrum with emission lines that can vary on a time- 
scale of 10 minutes jpuiimoto et alJllOOTh . while the diffuse ther- 
mal emission and CXB are expected to be more stable and easily 
modeled. At the low Galactic latitude of this source (6 = 3°), the 
soft thermal background could vary significantly with position. We 
opt to exclude the high energy end of the spectrum where the NXB 
component dominates with several emission lines, and to also ig- 
nore the low energy end where charge exchange and local thermal 
emission muddle the data, leaving the 1 — 5 keV band for consid- 
eration. Since the central regions have temperatures higher than the 
upper limit of the energy band used, we test that they are consis- 
tent with temperatures measured in a 1 — 8 keV band. Systematic 
uncertainties due to particle background subtraction at high ener- 
gies and Galactic and other low energy background components 
are discussed in Section |43] We opt for the narrower energy range 
to optimize the signal to noise for whole cluster. 

Since cluster emission fills the observed field, we cannot use 
detector regions from these pointings to subtract the remaining 
background. Instead, we analyze Suzaku data from the Lockman 
Hole (observation ID 102018010) taken only 8 days prior to the 
start of our observations. Despite the large difference in absorbing 
column density measured from H I maps (A^h = 5.7 x 10^® cm^'^ 
for the Lo ckman Hole versus 4.2 x 10^^ cm ~^ for PKS 0745-191 ; 
iKalberla e t al. 2005), ROSAT observations l lSnowden et ai]|l997l) 
detect similar background levels from 1 — 2 keV at the two posi- 
tions. We find that the Lockman Hole data from 0.5 — 8keV is 
well-modeled by a power law of photon index 1.4, absorbed by 
the column density measured by the H i maps, added to an un- 
absorbed MEKAL th ermal plasma model jMeweetal.|[l98l.ll986l : 
iLiedahl et al] Il995h at 0.1 keV with solar metal licity using the 
relative abundances of I Anders & Grevessd ( Il989t) , allowing only 
the normalizations to vary. The best fitting normalizations in the 
1 — 5keV range are 7.26 x 10~* photons keV~^ cm^^ s^^ at 
1 keV for the CXB and 0.78 cm~® arcmin"'^ for the soft thermal 
component0 We do not find an improvement in the fit with another 
low- temperature component, as is often used (e.g., IVikhlinin et al.l 
l200S j). We note that in the 1 — 5 keV band used here, the fit is 
not even particularly sensitive to the normalization of the remain- 
ing soft thermal component, but we include it as an empirical de- 
scription of the background model to be applied to the PKS 0745- 
191 data. We attempted to fit the outermost regions of the cluster 
observations with a similar model and an added low-temperature 
(~ 0.25 keV) Galactic component. The normalizations required 
are much higher than for the Lockman Hole and would be incon- 
sistent with the ROSAT observations for this source. The fit to the 

1 We follow ISato et alj <2007h in dividing the surface brightness nor- 
malization of the thermal component by the solid angle used in the 
ancillary response file generation, Q = 20', so that Norm = 
lQ-'20 J rienjjdV /{4:T7{1 + 2)^_D^}/f! in cm~^ arcmin"^, where 
Da is the angular diameter distance to the cluster. 



outer regions of PKS 0745-191 is improved with a higher tempera- 
ture component (> 1.5 keV), indicative of cluster emission cover- 
ing the entire field. 



3 SPECTRAL ANALYSIS 

Spectra were extracted in annular regions of radii 0' — 2.5', 2.5' — 
6', 6' - 9.5', 9.5' - 13.5', 13.5' - 18.5', and > 18.5', with the 
outermost region reaching to nearly 24'. For each annulus we de- 
fine an effective radius that is approximately an emission-weighted 
mean (McLaughlin 1993), r = [0.5(rf/^ + rlit)f^^. These re- 
gions are centered at right ascension and declination 07:47:31.325, 
-19:17:39.95 (J2000), which is coincident with the cD galaxy and 
both the peak and centroid of X-ra y emission. D ue to the low count 
rates in some channels, we use the lCashl (ll979') C statistic. 

We used XSPEC vI2.5.0 to simultaneously fit models of the 
spectra, with data from each pointing grouped by annulus unless 
specified otherwise. We first subtract the NXB spectra from iden- 
tical regions on the detector, scaled by the relative exposure time 
of the night Earth observations. Next, we model the local thermal 
emission and CXB components, fixed to the best-fitting values from 
the Lockman Hole observations. Finally, we model the remaining 
emission, that from the cluster, with another MEKAL component at 
the cluster's redshift, z — 0.1028. We fix the column density ab- 
sorbing the cluster and CXB components to the best-fitting value 
at the center, A'^h = 3.6 x 10^^ cm"'*, which is lower than the 
Hi measurement but consistent with the value from XMM's EPN 
jChen et al.l2003l) . We allow the metallicity to vary in the center but 
fix it to 0.3 Z© at larger radii. Of the remaining parameters, only the 
temperature and normalization of the thermal ICM component are 
allowed to vary. 

The broad point-spread function (PSF) of Suzaku, with a half- 
power diameter of ~ 2', will cause some emission from each annu- 
lus on the sky to be distributed into others on the detector. W e cor- 
rect for this effect following the method o f ISato et al.l f2007h . who 
importantly show that Suzaku's PSF is nearly energy-independent. 
We use a circularly sjonmetric /3-model with the best-fitting param- 
eters of Chen et £d] j2003l) . consistent with the surface brightness 
profile observed here, and a monochromatic energy of 1.5 keV to 
generate a photon list that is propagated with a ray-tracing simu- 
lator (xissim) through the detector optics. For each annulus on the 
sky, we use the ratio of flux landing in the corresponding detector 
annulus to the flux in each other annulus to calculate cross-annuli 
normalization factors. In the spectral model, we add a thermal com- 
ponent from each annulus to every other annulus with the temper- 
ature tied to the original region and the normalization scaled by 
these factors. The relative contributions to each annulus from other 
annuli, as derived from these PSF correction factors and the best 
fitting normalizations, are given in Table (2] The model was fitted 
simultaneously over all annuli and resulted in an excellent fit; the C 
statistic, though not a direct measure of goodness of fit, was 65413 
using 98640 bins and 98627 degrees of freedom. The spectra and 
models are shown in Figs. |2] and [S] We note that the individual 
spectra are fitted simultaneously, and are only combined by annu- 
lar regions in these figures for visual clarity. 
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0'-2.5' 


2.5'-6' 


6'-9.5' 


9.5'-13.5' 


13.5'-18.5' 


18.5'-24' 


0'-2.5' 


89.9 


19.3 


2.6 


0.4 


< 0.05 


0.1 


2.5'-6' 


9.7 


66.8 


5.7 


0.6 


0.1 


0.1 


6'-9.5' 


0.3 


13.1 


84.4 


7.8 


0.6 


0.1 


9.5'-13.5' 


0.1 


0.5 


6.5 


78.5 


7.3 


0.8 


13.5'-18.5' 


< 0.05 


0.1 


0.5 


11.9 


84.1 


10.9 


18.5'-24' 


< 0.05 


0.1 


0.3 


0.7 


7.8 


88.1 



Table 2. Relative flux contributions due to PSF spreading. Values con'espond to tlie fraction of flux measured in the detector region (columns) originating from 
the cluster region (rows), i.e., of the flux measured in the outermost annulus, we expect 88.1% to come from the corresponding region in the sky, and only 
0.1% to have originated in the central projected region. 




Energy (keV) 

Figure 2. Data and spectral fits including PSF-correction for each annular 
region, illustrating the softening of the spectra with increasing radius. The 
regions plotted are 0' - 2.5' (black), 2.5' - 6' (red), 6' - 9.5' (green), 
9.5' - 13.5' (blue), 13.5' - 18.5' (cyan), and > 18.5' (magenta). Data 
have been rebinned to a minimum significance of 5(t for each point, and 
each spectrum has been arbitrarily renormalized for visual clarity. 



4 RESULTS 

4.1 Temperature, density, and entropy profiles 

Temperature and density profiles are the primary data products of 
X-ray observations of clusters, from which other properties includ- 
ing entropy, mass, and gas fraction profiles can be derived. We 
obtain the density from the normalization of the best-fitting ther- 
mal plasma model which is proportional to the emission measure, 
J noTiH dV, where ric ~ 1.27ih for typical abundances. We use 
the simpifying assumption of constant density within each annu- 
lus, and estimate the volume as that of the intersection between 
cylindical and spherical shells with shared inner and outer radii, 
V — (47r/3)(rout — J'fn)'^^^, scaled by the fraction of the projected 
annulus observed. With the density and temperature determined, 
we can calculate the entropy profile, K = kT jri^J'^ . Each of these 
profiles is shown in Fig.|4] 

The cool core, though blurred by the wide PSF, is seen as a 
decrease in temperature near the center. The more novel result is 
the clearly-observed decline in temperature in the outskirts of the 
cluster. Excluding the core, this decline can be fit with a power law 
r(r) = Tor°^ with index qt = -0.94 ± 0.06. The gas density 



profile can similarly be fit with a power law, excluding the central 
region again for uniformity, with index a„ = —2.23 ± 0.01. If we 
assume a polytropic distribution for an ideal gas, P cx nnkT oc 
Wh, then we find T = 1 + {aT/a„) = 1.42 ± 0.03. This value is 
higher than that found for cooling flow clusters determined by e.g., 
De Gr andi & Molendi (2002) , who measured F ~ 1.2 in the radial 
range 0.2ri8o < r < O.Grigo. 

Gas with high entropy rises, so the radial entropy profile is ex- 
pected to increase. Puzzlingly, the entropy profile we observe rises 
from the center and then appears to level off beyond '-^ 10', putting 
it on the border of convective instability. We will discuss possible 
explanations for this behavior in later sections. 



4.2 Mass profiles 

To derive a mass profile, we assume that the gas properties are 
spherically symmetric and that the cluster is in hydrostatic equi- 
librium. Balancing the thermal pressure gradient against the gravi- 
tational potential, we obtain an expression for the total gravitating 
mass in terms of temperature and density: 



M(< r) 



dr 



dr 



(1) 



We assume that the total density distribution can be described by 
an NFW profile, 



p{r) 



Po 



r/rs (1 + (r/rs)) 



(2) 



where po is a density normalization and is a scale radius related 

to the cluster's concentration, C200 = f'20o/fs. 

Following an approach similar to that of ISchmidt & AUetil 

we use an NFW model and the observed gas density dis- 
tribution to predict a temperature in each annulus, beginning with 
the outermost bin. Density gradients are taken linearly between the 
radial bins, which reduces the correlation in uncertainties between 
non-adjacent data points that are created with oth er interpolations 
such as cubic splines, as some authors use (see IVoigt & FabianI 
1 20061 for a discussion). This approach could be improved with 
more radial bins, but the wide PSF of Suzaku and the low count 
rate in the cluster outskirts limit us from using thinner annular re- 
gions. To ensure that the outermost temperature value is robustly 
estimated, we take the median outer value of the profiles that best 
fit 100 Monte Carlo realizations of the observed temperature pro- 
file, normally distributed within its error bars. The model tempera- 
ture profiles are similar if we begin at the innermost region, so the 
result is not particularly sensitive to the boundary values. Iterating 
over a range of NFW parameters, we select the mass model that 
produces the temperatures that best fit the observed profile. 
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Figure 3. Same as Fig. |2] but with residuals to show the quality of the 
spectral fits, and with the Lockman Hole background spectrum included for 
comparison. Count levels have not been renormalized and we note that the 
regions have different areas and exposure times. Dotted and dashed curves 
represent the background (CXB+thermal) and cluster emission models, re- 
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Figure 4. Projected temperature, density, and entropy profiles. The vertical 
dashed line shows our estimate of r200, derived in Section l4!2l 



In order to resolve the region within the NFW scale radius, we 
supplement the Suzaku data with an archival 20 ks Chandra obser- 
vation (ID 2427) taken with the ACIS-S detector in VF AINT mode 
These data were reduced following the description of lO'Dea et alj 
(2008'), with spectra extracted in circular annuli and deprojected 
using the method of Sanders & Fabian (2007). After using the out- 
ermost Chandra annulus to subtract contaminating flux from in- 
ner regions it is omitted from further analysis since it may contain 
external cluster emission, leaving six regions to replace the inner- 
most Suzaku bin, where the temperature and density is in reason- 
able agreement with the Chandra data. We use C handra data for 
its sha rper PSF, but the XMM temperature profile in lSnowden et al.l 
( l2008l) is also consistent at the smaller radii of that observation. 

The model temperature profile, shown in Fig. |5] does not fit 
the data very well over all annuli, and the best fit to the observed 
temperatures (x^ ~ 39 for 7 degrees of freedom) is obtained after 
ignoring the outermost Suzaku bin which is beyond the virial ra- 
dius where the assumption of hydrostatic equilibirium is expected 
to fail, and the Chandra regions in the cool core within ~ 70 kpc of 
the center where the gas is likely to be multiply-phased. The main 
deviation in our data from this acceptable temperature fit arises 
from the 2.5' — 6' Suzaku annulus. Possible explanations include 
a PSF-correction that could be too large, a cross-calibration issue 
between Chandra and Suzaku, or even non-gravitational heating in 
that region. 

The resulting mass profile has best-fitting NFW parameters of 
C200 ~ 7.5lj g and Ts — 230I4Q kpc, producing a virial radius 
''200 = C2oo»'s = 1.72 ± 0.06 Mpc, or 15.2' and mass A/200 = 
200pc(z)(47r/3)(r2oo)^ = 6.4 ± 0.6 x lO" M© . We estimate 
the total uncertainty on the virial mass, including systematics, of 
30%, which is described further in the next section. We compare 
these values with previous determinations in Table [3] and plot the 
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Figure 5. Chandra (red triangles) and Suzakii (black crosses) temperature 
values used in mass fitting analysis, along with matching temperatures for 
best fitting NFW profile (blue circles). 



Figure 7. Gas fraction profile for PKS 0745-191, the ratio of the dotted to 
solid curves in Fig.|6] with Icr statistical eiTor range in gray. Dotted lines 
show the Icr range around the value of Cm from Komatsu et al. ( 200^. 
The vertical dashed line shows our estimate of r200> derived in Section l431 
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Figure 6. Cumulative profiles of gas mass (dotted lower curve) and total 
gravitational mass (solid upper curve) with Icr statistical error ranges in 
gray. The single eiTor bar shows the 30% systematic uncertainty on Af2oo 
and the vertical dashed line shows our estimate of r2oo. derived in Sec- 
tionl42l 



mass profile in Fig. |6l along with the cumulative gas mass profile 
obtained from the power-law fit to the Suzaku gas density profile. 

The NFW concentration is higher than previous estimates, and 
the virial mass we determine i s lower, due in part to the sm aller 
virial radius found in this work. lReiprich & B6hringe3 j2002h . who 
used an isothermal fit to ROSAT and ASCA data and found a value 
of 1.95 ± 0.03 Mpc updated to the cosmological parameters we 



adopt. ISchmidt & AUer] tOOlh found an even larger value for the 
virial radius, fitting an NFW profile to Chandra data in the central 
region to obtain r2oo ~ 2.2 Mpc. One issue with past NFW fits 
used to calculate the virial radius is that the data often do not extend 
even to the NFW scale radius, let alone the virial radius. 

We can compare our mass profile with the projected mass de- 
termined from a str ong lensing arc of radius 34.5 kpc seen with 
Hubble imagery by lAllen et al.l Updating their mass esti- 

mate for the circular lens model to account for the current cosmo- 
logical parameters, we expect a mass of 6.4 x lO^'^ Mq projected 
within a cylinder along the line of sight of the center of the cluster 
with a radius of the lensing arc, which is consistent with the sta- 
tistical and systematic uncertainties of our best-fitting NFW model 
when integrating the mass in the cylinder out to r2oo- We note that 
this includes an extrapolation of our profile into the core region 
where we have excluded data from our fit because of nonthermal 
processes. While it is reassuring to see this reasonable agreement 
between separate mass estimates, we caution that comparisons of 
this type of wide-field X-ray analysis with strong lensing data only 
involve a small fraction of the total cluster mass. Weak lensing anal- 
yses, which can provide mass projections over a much broader area, 
will provide important comparisons for virial mass estimates. 

In Fig.|7l we plot the gas fraction profile, defined as the frac- 
tion of mass within a given radius composed of X-ray emitting gas. 
The gas fraction appears to rise slightly at the outer radii, consis- 
tent with the profi le observed for PKS 0745-191 in ROSAT data by 
lAllenetalJ ( ll996t) . The values of the /gas profile in the outskirts of 
this cluster are s imilar to recent meas urements of the mean cosmic 
baryon fraction jKomatsu et al.ll20o3) . 



C200 



M2oo(lO^*M0) Reference 



3.83 
5.12 
5.46 
5.86 
7.5+ 



+0.52 
0.27 
+0.40 
0.40 
+3.22 
-2.88 
+ 1.56 
-1.07 
1.1 
0.9 



+3.5 
4.0 
+ 1.2 
1.2 

+52.2 

+0.6 
0.6 



18.6 
10.0 
9.7 



6.4"' 



Allen et al. (2003) 
Pointecouteau et al. (2005) 
Voig t & Fabian (200&) 
Schmidt & Allen (2007) 
This work 



Table 3. Conce ntration and virial mass determ inations for PKS 0745- 
191, compiled bv lComerford & NataraianI j2007l) with values for this work 
added. 



4.3 Uncertainties 

A number of systematic uncertainties could add contributions to 
our error budget beyond the statistical ranges plotted. We sum- 
marize our estimates of several of these systematic uncertainties 
to the temperature profile, including the effects of PSF-correction 
and renormalization of background components, in Table |4] The 
changes due to the PSF-correction are mostly small, with the main 
effect being a 12 per cent increase in temperature in the second 
innermost annulus after accounting for emission from the cool gas 
in the core that has scattered outward. There is a significant un- 
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Radius 


kT 


Stat. 


NoPSF 


NXB 


CXB 


Gak 


TVh Z 


(') 


(keV) 






±3% 


±10% 


±50% 


±10% 0.2- 0.4 Z0 



0-2.5 


6.92 


+0.12 
-0.12 


+0.08 


< 0.005 


< 0.005 


< 0.005 


0.68 


< 0.005 


2.5-6.0 


7.85 


+0.45 
-0.39 


-0.82 


< 0.005 


< 0.005 


0.01 


0.90 


0.02 


6.0-9.5 


4.90 


+0.30 
-0.27 


-0.01 


0.02 


0.04 


0.03 


0.41 


0.05 


9.5-13.5 


3.22 


+0.47 
-0.36 


-0.05 


0.14 


0.14 


0.06 


0.28 


0.07 


13.5-18.5 


2.15 


+0.19 
-0.18 


+0.09 


0.06 


0.17 


0.06 


0.15 


0.11 


> 18.5 


2.07 


+0.25 
-0.23 


+0.02 


0.10 


0.21 


0.06 


0.13 


0.12 



Table 4. Estimated uncertainties in the radial temperature profile due to removing PSF-corrections or changing normalizations of background components and 
spectral model parameters within the ranges given. Each parameter is varied independently and we present only the maximum of the upward and downward 
temperature shifts (in keV) for clarity. Best-fitting PSF-corrected temperatures and their statistical uncertainties are included for comparison. 



certainty in the temperature of the central region be cause Suzaku's 
PSF blurs the steep gradient of the cool core. As Reiprich et al.l 
( l2008i) points out, the approach to PSF-corrections described in the 
previous section does not fully account for this blurring within the 
central region, however our method provides a reasonable approx- 
imation to the corrections needed. 

To ensure that the background models used do not signifi- 
cantly impact the results, we vary the normalizations of these com- 
ponents and measure the deviations produced in temperature. The 
uncertainty in the NXB is ~ 3 per cent (jTa wa et al. 2008), a nd from 
an analysis of ROSAT data by ICarrera. Fabian & Barconsl jl997h . 
we estimate that the cosmic variance due to point source clustering 
in the CXB is ~ 10 per cent over this field of view. We also test 
the effect of metallicity from 0.2 — 0.4 Z© in the outer annuli and 
a range of ±10 percent in the column density. Each of these vari- 
ations results in a change in temperature of < 10 per cent for all 
annuli, and the individual systematic uncertainties are smaller than 
the Poisson uncertainties for all but the column density's effect on 
the central regions where the count rate is high. Finally, the normal- 
ization of the local thermal component, which could vary the most 
across the cluster field of view because of the low Galactic latitude, 
actually has a very small effect on temperature fits above 1 keV. 

Without spectral deprojection, the data from each annulus 
will be contaminated by emission from regions at larger radii. 
But because the density profile falls so steeply and the ther- 
mal bremsstrahlung emission scales as the square of the density, 
the only noticeable eff ect is in the cluster core (e.g., Fig. 3 of 
ISanders & Fabiarj|2007h . This region is well within our innermost 
radial bin, so we do not expect significant errors due to the lack of 
deprojection of the Suzaku data. 

As an additional check on the reliability of our temperature 
and density profiles, we compare the values from Suzaku with the 
Chandra observations used in the mass estimate. The high spa- 
tial resolution of Chandra allows for our central region to be di- 
vided into several annuli, with the average temperature in agree- 
ment with the Suzaku value. As shown in Fig.jS] the density profiles 
are also reasonably consistent between the two observatories, with 
a slightly steepening slope at the larger radii seen by Suzaku. 

We can check the assumption of spherical symmetry used 
when deriving the mass profile, seeing if the emission is at least 
circularly symmetric by considering the projected profiles in each 
of the pointings separately. Fig. [9] shows the temperature, density, 
and entropy profiles for each pointing. Here we do not include PSF- 
corrections in order to keep each set of observations independent. 
The uncertainties are larger when splitting up the data, but the re- 
sults from the separate pointings broadly overlap, and no single 
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Figure 8. Comparison between Chandra (red triangles) and Suzaku (black 
crosses) density profiles. Error bars are smaller than the plot symbols. 



direction is consistently above or below the average temperature 
or density for every annulus. Scatter that is larger than the statisti- 
cal uncertainties could be attributed to remaining point sources or 
minor asymmetry in the cluster. 

The entropy seen in the outermost radial bin of the NW point- 
ing is significantly lower than nearby regions, mainly due to the 
temperature drop seen there. However, scatter between pointings 
or individual outliers are unable to explain the observed flattening 
of the entropy profile in the summed annuli. We plot a power law 
rising as r^'^ to compare with the entropy profile predicted by an- 
alytical models of accretion shock heating and seen in some obser- 
vations as well as numerical simulations jTozzi & Normanll200ll: 
IPonman et a l.'2003VVoit et al.'2005h. A smaU number of objects in 
the sample of Ponman et al. ( 2003) have entropy curves that do not 
rise with radius. The observed and expected profiles agree at small 
radii, rising out of the core, but there is a clear entropy deficit in the 
outskirts of the cluster, perhaps indicative of infalling gas which is 
not dynamically stable. Our background model would have to be 
off by an order of magnitude, or the column density changed by a 
factor of 2, to increase the entropy in the outskirts to the predicted 
values. 

For the estimates of the virial radius and mass, the method 
of fitting an NFW model to our temperature profile is an additional 
source of uncertainty. We have excluded the inner and outer regions 
where the assumption of hydrostatic equilibrium is likely to break 
down, but several variations of the temperature profile with these 
points included or excluded produce values for r20D consistent to 
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Figure 9. Comparison of temperature, density, and entropy profiles in each 
direction. The temperature in the NE direction is not well-constrained from 
9.5' — 13.5' due to a large amount of area excluded because of point 
sources, so the temperature and entropy points are not plotted in this re- 
gion. Slight offsets are added to bin radii for clarity. The dotted curve in the 
bottom panel shows K oc ' ^ . The veilical dashed line shows our estimate 
of ?'200: derived in Section l4!2l 



~ 10 per cent. Similar values are also obtained when using only 
the Suzaku data, which has a large uncertainty in the central tem- 
perature due to PSF-blurring. The wide radial span of this data set 
provides significantly tighter constraints on the NFW parameters 
than could be obtained with the Chandra data alone, as shown in 
Fig. [To] Data from the cluster core are still excluded in the con- 
tours plotted, though they might be used in a typical analysis. In 
either case with only Chandra data, the NFW scale radius is not 
well-constrained since the data do not extend out that far, and the 
resulting estimate for r2oo, and thus A/200, would be significantly 
larger than the values we obtain. We can also include the systematic 
uncertainties in the temperature profile within this process of fitting 
over NFW parameter space. Even with differences of 0.5 — 1 keV 
between the FI and BI sensors which have similarly shaped but off- 
set temperature profiles, changes in the best fitting value of i?20o 
are of order 10%, giving a systematic uncertainty of 30% for the 
determination of A/2oo • 

We can also estimate the cluster mass independently of any 
model parametrization by plugging our density and temperature 
profiles directly into Equation[T]and calculating gradients from the 
finite difference quotients of adjacent points. At large radii, the cu- 
mulative mass estimate using this method actually decreases, indi- 
cating that the hydrostatic assumption fails, and the uncertainties 
are large because spatial gradients are taken between points sepa- 
rated by the wide thickness of annular bins. However, at the virial 
radius near 15', the mass estimate agrees with that derived from the 
NFW model to within ~ 30 per cent. 




Concentration 

Figure 10. Likelihood contours for the concentration and scale radius pa- 
rameters of the NFW fits, drawn at Ax^ = 2.30, 6.17, and 11.80, corre- 
sponding to the Icr, 2(7, and 3(t confidence intervals, respectively. Thick 
red contours show the constraints from the combination of Chandra and 
Suzaku data, excluding the cluster core and the outermost annulus which is 
beyond the virial radius. Thin blue contours show the constraints without 
the Suzaku data for comparison. 



5 DISCUSSION 

These are the first well-constrained measurements of the gas tem- 
perature profile beyond r2oo for a rich and relaxed cluster, which 
we can now compare with model predictions and simulations. We 
see that there is a clear decline in temperature between 0.3 — 
1 . 5r-2 oo, consistent with recent observational trends at growini 
radii jVikhlinin et"ai]|2005L [Pratt et al.ll2007l : iReiprich et alj|200l 
and suggesting that conduc tion does not pla y a sig nificant role on 
typical cluster time-scales. iRoncarelU etal] (12006!) measured the 
magnitude of the temperature drop in hydrodynamical simulations, 
finding a 40 per cent decline f rom 0.3 — lr2oo, the same size found 
in the analytical treatment of lOstriker. Bode & Babull ( l2005i) . The 
drop observed in PKS 0745-191 over this radial range is closer to 
70 per cent, but the general shapes of the profiles are similar. 

The observe d entropy profile deviate s from expectations near 
the virial radius. iTozzi & Noiiiiaril ( 1200 ih describe the effects on 
entropy of the accretion phases, including adiabatic infall, shock 
heating, and adiabatic compression followed by cooling within the 
halo. Interestingly, we see evidence for an accretion shock beyond 
the virial radius in one direction (NW pointing in Fig. |9), sugges- 
tive of cool material falling inward along a filament. The change in 
entropy in the outer annulus of this region is not seen in the other 
directions, though the values there are still lower than would be 
expected from the typical increase of K oc r^'^. 

The sharper decline in temperature than predicted and the flat- 
tening of the entropy profile in the outskirts suggest a need for 
nonthermal pressure support in order to maintain dynamic stability. 
The numerical simulations of Eke, Navarro & Frenk ( 1998) show a 
significant rise in the ratio of bulk kinetic energy to thermal energy 
beginning within the virial radius. Merger activity increases the tur- 
bulent pressure, and even though PKS 0745-191 appears relaxed 
morphologically, it should still be accreting matter in the outskirts. 
We note that the infall time-scale, tt ~ {r^ / GMY^'^ , is of order 
1 Gyr at the virial radius, and the time-scale set by the speed of 
sou nd, ts (r ^ mn/r fcT)^''^, is several times larger. 

iNeumannl ( |2005|) . studying the summed profiles of 14 nearby 
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Abell clusters beyond r2oo, argues that the ICM at large radii 
may not be in hydrostatic equilibrium, and that cool gas not seen 
in X-rays could increasingly d ominate the baryon content in the 
outskirts. lAfshordi et alj |2007) stacked WMAP observations of a 
large sample of massive clusters and found a deficit in thermal en- 
ergy in the outskirts from the Sunyaev-Zel'dovich profile, also ar- 
guing for cool phase of the ICM. It is not well understood how these 
different phases would mix, and complicated gas physics would 
likely result. The assumption of hydrostatic equilibrium fails be- 
yond the virial radius, as seen by the poorer temperature fits from 
mass modeling when the outermost temperature value is included. 
It is not obvious to what extent the ICM within the virial radius is 
not in equilibrium, but the entropy and temperature profiles sug- 
gest that nonthermal physics may be important in the outskirts of 
this cluster. 

In Fig.[TT] we plot the observed temperature profile against 
two models. The first is the expected temperature if the gas content 
traced the NFW distribution of dark matter, with a Keplerian veloc- 
ity profile, i.e., T{r) oc v{r)^ oc M(< r)/r. The second and more 
realistic case is for the projected X-ray temperature of gas arranged 
polytropically and in hydrosta tic equilibrium with an N FW gravita- 
tional potential, as derived bv iSuto. Sasaki & Makind {l998). The 
outer temperatures we have observed are lower than those expected 
for reasonable values of the model parameter^- Additionally, we 
can fit the observed temperature profile using a simple parametriza- 
tion 

r(r)=ri(^^^ (l + f)''keV (3) 

with best-fitting parameters Ti — 135, a = 2.4, f3 = —4.1, where 
rs is the scale radius from the best-fitting NFW profile. We note 
that the normalization is related to the value of the temperature at 
this radius, Ti = T{r^)/2'^. 

A possible explanation for the deficit of thermal energy seen 
at large radii is that infalling matter has retained some of its kinetic 
energy in bulk motion. If there is increasing pressure support from 
bulk motion, turbulence, or other nonthermal processes in the outer 
regions of the cluster, then the gravitational potential and thus total 
mass would be incorrectly estimated when assuming hydrostatic 
equilibrium. This effect would depend on the changes in both the 
temperature and its gradient, making estimation of the true mass 
difficult. 

Alternatively, the peak temperature, seen in the second- 
innermost armulus, could be inflated due to nonthermal processes 
such as shock heating during infall. AGN heating is thought to 
play an important role in cool core clusters, and some excess en- 
ergy could also heat the gas beyond the core. More observations 
are needed of this cluster and others to determine how well the 
self-similar scaling relations, which depend on the dominance of 
gravitational process in creating equilibrium, apply out to the virial 
radius. 

The gas fraction we obtain is more consistent with previous 
results. At r2oo, /gas = 0.13, lower than that found with XMM 

^ We set the parameter Bp, defined by ISuto et al J )l998h . equal to unity, 
and use the polytropic index derived from the best-fitting power laws to our 
data excluding the central region, T = 1.5. It is possible to obtain tem- 
perature profiles that decline as steeply as the observed one by increasing 
the value of Bp, though this would require tmncating the gas distribution 
shortly beyond the regions we have measured. For the range F = 1.1 — 2, 
the polytropic index used does not have a significant effect on the shape of 
the scaled profile. 
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Figure 11. Suzaku temperature profile, normalized by the emission- 
weighted mean and compai'ed with model profiles. The solid curve shows 
the parametrization of Equation [5] also rescaled by the emission- weighted 
mean temperature that we observe. The dotted curve shows the tempera- 
ture for gas tracing an NFW distribution, while the dashed curve shows a 
polytropic profile in hydrostatic equilibrium with an NFW potential. Both 
models are normalized to match the maximum observed temperature. 



IChe n et al.ll2003h . but similar to the ROSAT value of lAUen et all 
( 19961). The trend of /gas rising with radius is also seen in other 
clusters and in simulations (e.g., IVikhlinin et al.ll2006l : Iske et al.l 
IT998 ), since the gas density distribution is less centrally concen- 
trated than that for collisionless dark matter. The outermost mea- 
sured value of /gas is consistent with the mean cosmic baryon ra- 
tio derived from WMAP da ta and other cosmolo gical observations, 
Slb/f^m = 0.165 ± 0.006 dKomatsu et alj2008l) . 



6 SUMMARY 

We have presented the first measurements of the gas beyond the 
virial radius from a rich and apparently relaxed galaxy cluster. The 
temperature profile shows a significant decline at outer radii, while 
the entropy profile levels off. The spatially resolved temperature 
and density measurements out to and beyond the virial radius pro- 
duce improved constraints on the mass and gas fraction profiles, 
though the observed temperature profile does show deviations from 
that expected from cluster gas in hydrostatic equilibrium with an 
NFW mass profile. 

Observations of the outskirts of clusters offer a direct probe 
of their assembly history. The apparent shock front in one direction 
makes it clear that this cluster is still accreting material, likely along 
a filament. Evidence for nonthermal pressure support suggests that 
bulk motions from merger activity could be making a significant 
contribution to the gas energy in the outskirts of this cluster. 

PKS 0745-191 is an ideal candidate for such observations, 
given its distance and luminosity. The angular size fits reasonably 
within a few pointings and the count rate is high enough for good 
spectral statistics. The low and well-constrained background lev- 
els of Suzaku are crucial for low surface brightness m easurements 
like those presented here and bv lReiprich et alj ( l2008h . We expect 
that as the number of cluster observations in the Suzaku archive 
increases, these types of analyses will be able to probe the self- 
similarity of the outskirts of clusters. More studies at the virial ra- 
dius, including X-ray observations, gravitational lensing, the SZ 
effect, and numerical simulations, will provide a better understand- 



© 0000 RAS, MNRAS 000, 000-000 



10 M. R. George et al. 



ing of the regime where hydrostatic equilibrium breaks down and 
clusters are still being assembled. 
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